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Turbulence Effects on Optical Beams f A6A 


• Scintillation 

• Beam broadening 

• Spatial coherence 

• Angle of arrival 

• Temporal pulse 
stretching 


Distorted 



Temporal and spatial intensity fluctuations at the receiving 
aperture results in power surges and fades 
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Turbulence Characterization r** 



Requires simultaneous 
measurements at f and f + A f 


Temperature (or humidity) structure 
function 
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Contains the spatial statistics of the 
temperature field 


Within a range of certain A r, the 
well-known Kolmogorov “2/3” law 
holds , 
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Single Radiometer Turbulence 
Characterization 



Wind velocity field 
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1 st measurement 


2 nd measurement 


• Temperature (or humidity) 
measurements are taken at a 
fixed location 

• Wind velocity field shifts the 
turbulent air mass by a 

distance A r = UAt between 
measurements 


• For a particular altitude, the 
structure function is now 


D T (t,t + UAt) = / 

v 

T(t) - T(t + UAt) 

’> 



Due to the radiometer integration 
time At, Kolmogorov-Obukhov 
turbulence theory and the Taylor 
frozen flow hypothesis must be 
modified 
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Energy Transfer Spectra 



• From the Boussinesq approximation (eddy viscosity model) of the Navier- 


Stokes equations, it is possible to obtain equations involving the Fourier 
spectra of the turbulent energy, wind velocity and temperature 
fluctuations* 


Energy spectrum 

related to Reynolds Buoyancy 
stresses 


Kinematic 

viscosity 


Turbulent energy 
fluctuations spectrum 


Energy transfer 
spectra due to 
wind velocity 


F(k) - 4>1 3 (fc)f + / 34> 3T (k ) - 2 vk 2 4>{k) = 0 
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Energy transfer spectra 
due to temperature 
fluctuations 


Energy spectrum 
related to heat flux 


Thermal 

diffusivity 


Temperature 
fluctuations spectrum 


)|C 

Tchen, C.M., "On the Spectrum of Energy in Turbulent Shear Flow," J. Res. Natl. Bureau of Standards 50 (1), 51-62 (1953) 
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Modified Turbulence Spectrum 


Case 1: Near the boundary 
surface 

- Significant stratification and shear 
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Case 2: Free atmosphere 

- No stratification or shear 


Heisenberg 
constant = 1 
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TKE 

dissipation 


Thermal diffusivity/ 
kinematic viscosity 


Thermal 

dissipation 


(fyrri.k') = Bk 

B = Q) 4 1 /3 7 - 2 /3 6 -l A f e -l / 3 


General: 

- Asymptotically reduces to 
either Case 1 or Case 2 
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Relating to the Structure Functions 



*00 


—00 


D t ( A t) = 2 ( 1 - ( e~ ik< - u+v l At ) ) V TT (k)dk 


Average wind 
velocity 


Wind velocity 
fluctuations 


The connection of the 
temporal statistics of the 
temperature T(t) to the 
spatial spectrum VttW 
is through the Fourier- 
Stieltjes transform 


Evaluation of the 
integral is analytical in 
terms of Meijer G 
functions, however 
two useful series 
expansions can be 
obtained for the 
asymptotic cases 
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Kolmogorov Crossover frequency 

"2/3" law kc = (B/A) m 
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Instrumentation 




NASA TDRSS ground terminal site 
located at White Sands, NM 

Radiometrics MP-3000A 

35 calibrated channels 

- 300 MHz bandwidth/channel 

- 21 K-band (22 to 30 GHz) 

- 14 V-band (51 to 59 GHz) 

1.1 second integration time per 
channel 

Total At ~ 40 second sample 
period 

Temperature resolution ~ 0.1 K 
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Temperature (K) 


Temperature Data Analysis 


Example fluctuations over a 2-hour time period 



Fluctuation 
standard deviation, 
a = 0.5 K- 3 K 


Average profiles for each 10 min window 



Temperature (K) 


Average lapse rate 
of = -5.474 K/km 


Measurements taken 
in January, 2013 

Dataset comprised of 
about 2100 profiles 
taken over a 24-hour 
period 

Each temperature time 
series divided into 10- 
minute moving average 
windows 
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Statistical Wind Model 



Vertical profiles of horizontal 
wind speed 

SPARC Data Center High- 
resolution radiosonde 
measurements at Santa 
Theresa, NM 

Statistics derived from 2376 
wind profiles 

Principal component analysis 
(PCA) used for data reduction 
and retention of key features 
of the wind behavior 



Ground wind speeds 
typically between 5 
and 8 m/s 


Max tropopause 
wind speeds 
typically between 
30 and 40 m/s 


Statistical Wind Model Profiles 
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Results for C n 2 



Refractive index structure parameter profiles 



10 15 10' 14 10' 13 
Structure factor, C „ 2 



General Hufnagel-Valley model 


• At optical wavelengths the 
refractive index structure 
parameter is a function of only 


r 
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Average prevailing 
pressure (hPa) 


Free atmospheric 
estimates of C„ about 
lOx - lOOx larger than 
expected 


Ground estimates 
of Cl = 10 13 m' 2/3 


Specific atmospheric 
conditions during 
data compilation 
were not available, 
thus a nominal value 
of k c = 15 m -1 was 
assumed 
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Results - Coherence Diameter 



Coherence Diameter vs Time 
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Night time Day time variation 

variation between between 4.3 cm - 6 cm 
5.2 cm - 7.5 cm 


• Coherence diameter, also known 
as the Fried parameter 



• Determines resolution limitations 
of telescopes 


D > r 0 

atmosphere limited 


D < r 0 

diffraction limited 


• Also determines the spacing of 
actuators in adaptive optical 
systems 
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Results - Greenwood Frequency 



The Greenwood frequency 
specifies the response 
characteristic required of an 
atmospheric adaptive optics 
system to mitigate the refractive 
index perturbations 




These values are about a factor of 
3 larger than expected for the 
experimental site in January 


Night time Day time variation 

variation between between 200 - 350 Hz 
150 - 290 Hz 
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Results - Five Days in October 2012 


Diurnal variations are 
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Summary 



• Atmospheric remote sensing method using a single microwave 
profiling radiometer to obtain temperature and humidity 
turbulence structure parameters 

• Augmented Kolmogorov turbulence theory to account for boundary 
effects in a general stratified atmosphere 

• Test case shows promising results; however Greenwood 
frequencies and coherence diameters are over/under estimated 

• A more rigorous turbulence spectrum derivation is required 

• Ground-based measurements of the gradient Richardson numbers 
are required for better estimation of the crossover frequency k c 

• Concurrent radiosonde measurements of structure parameters 
along with the radiometer is needed for appropriate comparison 
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